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ABSTRACT 


In  an  earlier  report,  a new  class  of  multilayer  bounded 
cellular  automata  was  defined  and  shown  to  improve  the  lower 
bound  recognition  time  for  many  basic  array  recognition  tasks. 
We  continue  our  investigation  of  pyramid  cellular  acceptors 
by  presenting  new  results  on  their  capabilities  as  two- 
dimensional  pattern-recognizing  machines. 
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1. 


Introduction 


A central  concern  of  image  analysis  is  to  construct 
algorithms  for  extracting  a description  from  a given  picture. 

In  general,  this  involves  segmentation  of  the  picture  into 
parts,  measurement  of  properties  of  the  parts,  and  determina- 
tion of  relations  among  the  parts.  To  this  end,  many  basic 
picture  analysis  techniques  have  been  developed.  For  example, 
thresholding,  template  matching,  edge  detection  and  connected 
component  analysis  are  a few  of  these  techniques. 

The  design  and  evaluation  of  efficient  image  analysis 
techniques  and  algorithms  have  not  been  adequately  researched. 
In  particular,  the  tradeoffs  between  computer  architecture  and 
algorithm  complexity  for  a given  problem  are  of  considerable 
interest  because  of  the  high  data  rates  necessary  in  many 
image  application  areas,  while  hardware  costs  continue  to  drop. 
Furthermore,  many  basic  image  analysis  operations  operate 
independently  on  each  point  of  the  image  and  its  neighbors. 
Cellular  parallelism  is  very  efficient  in  computing  such  local 
operations.  In  fact,  the  most  obvious  technique  computes  each 
point's  new  value  independently,  leading  to  a very  inefficient 
serial  algorithm  (due  to  neighborhood  overlap) , while  being 
ideally  suited  for  implementation  as  a parallel  algorithm. 

In  parallel  computations  a primary  concern  is  the  distri- 
bution of  results  aunong  the  processors.  The  arrangement  of  the 
data  and  the  processors  handling  them  must  be  considered  in 
light  of  complicated  data  flow  that  may  be  required.  This 


problem  is  especially  evident  in  cellular  architectures,  in 
which  processors  are  uniformly  interconnected  to  a fixed 
number  of  neighbor  processors.  Under  conditions  where  communi- 
cation between  distant  processors  is  necessary,  unavoidable 
inefficiencies  may  arise. 

With  these  factors  in  mind,  we  introduced  in  [1]  a pyramid 
cellular  machine,  in  which  the  processors,  or  cells,  are  con- 
figured in  a logarithmically  tapering  pyreunid.  Each  pixel 
in  the  input  picture  is  associated  with  its  own  cell  in  the 
base  of  the  machine.  Each  cell  operates  on  its  own  local 
memory  and  has  its  own  local  control,  which  is  globally  syn- 
chronized by  a discrete-step  clock.  Each  cell  can  examine 
the  memory  contents  of  its  nine  neighbors  — its  "father"  cell 
in  the  layer  above,  its  four  nearest  neighbors  in  its  own  layer, 
and  its  four  "son"  cells  in  a two-by-two  block  in  the  layer 
below.  The  model  is  defined  by  stacking  cellular  array  auto- 
mata, where  each  layer  contains  one  quarter  the  number  of 
in  the  layer  below  it. 


processors 


* 

In  [1]  this  new  model's  computing  power  was  compared 
with  that  of  several  others,  and  some  basic  pyraunid  recognition 
algorithms  were  described.  It  was  shown  that  pyreunid  cellular 
acceptors  can  solve  many  nontrivial  tasks  in  time  proportional 
to  the  logarithm  of  the  diameter  of  the  input.  This  is  in 
constrast  to  conventional  bounded  cellular  acceptors  which 
require  at  least  time  proportional  to  the  diameter  of  the 
input  for  recognition.  For  the  most  part,  however,  the 
recognition  results  presented  in  [1]  were  confined  to  the  one 
dimensional  case,  i.e.,  when  the  pyramid  reduces  to  a tri- 
angle of  processors  and  the  input  is  a string  of  symbols. 

This  report  extends  the  analysis  of  cellular  pyramids  by 
providing  new  results  on  recognizing  two-dimensional  picture 
languages.  We  first  review  the  definition  of  a cellular 
pyramid  in  Section  2.  Section  3 describes  some  basic  pyramid 
recognition  algorithms,  including  local  property  detection  and 
counting,  palindromes,  rectangles  and  squares.  In  Section  4 
we  compare  pyramid  acceptors  with  finite-state  acceptors  and 
bounded  cellular  acceptors. 

« 

; The  cellular  pyramids  defined  in  this  report  are,  in 

general,  nondeterministic ; but  the  algorithms  given  in  this 
;•  report  are  all  deterministic.  Some  results  on  nondeterministic 

r cellular  pyramids  will  be  presented  in  a subsequent  report. 


2. 


Definitions 


In  this  section  we  review  the  definitions  of  a bounded 
cellular  array  acceptor,  pyramid  cellular  acceptor,  and 
bottom-up  pyramid  cellular  acceptor.  We  also  introduce,  for 
comparative  purposes  in  Section  4,  the  up-down  pyramid  cellular 
acceptor. 

A bounded  cellular  array  acceptor  (CA)  is  a finite, 

rectangular  array  of  identical  finite  state  machines  (FSM's), 

or  cells.  Each  of  these  cells  is  a quadruple  M=  (Qj^,Q^,  6 ,A)  , 

where  is  a nonempty,  finite  set  of  states,  c is  a 

nonempty,  finite  set  of  input  states,  and  A t is  a nonempty, 

5 % 

finite  set  of  accept,  or  final,  states.  2 is  the 

state  transition  function,  mapping  the  current  states  of  M 
and  its  four  nearest  neighbors  into  a set  of  possible  next 
states  for  M.  If  M is  deterministic,  the  range  of  6 is  a 
single  state  from  Associated  with  the  CA  is  a special 

boundary  state  Qjj  which  can  never  be  created  or  destroyed 
by  6 . 


A configuration  of  a CA  is  an  assignment  of  states  from 
to  each  cell  in  the  CA.  A step  of  computation  consists 
of  a single  application  of  the  state  transition  function  5 
simultaneously  at  each  cell.  An  input  configuration  is  a 
configuration  before  the  first  step.  If  a cell  is  on  the 
border  of  the  array  its  initial  state  must  be  the  boundary 
state  #,  otherwise  its  state  may  be  chosen  from  Q^.  Cells 
initially  in  the  boundary  state  are  called  boundary  cells; 


all  others  are  called  retina  cells.  An  input  configuration 
is  accepted  by  a CA  if  at  some  step  the  upper-left  corner 
retina  cell  enters  an  accept  state.  The  set  of  input  con- 
figurations accepted  by  a given  CA  defines  its  language, 

A pyramid  cellular  acceptor  (PA)  is  a pyramidal  stack 
of  CA's,  where  if  the  bottom  array  has  retina  size  2 by  2 , 

— 1 IT*  1 

then  the  next  lowest  array  has  size  2 by  2 , and  so  on, 

until  the  (r+l)st  layer  consists  of  a 1 by  1 CA.  This  apex 
cell  in  a PA  is  called  the  root.  Each  cell  is  an  identical 
FSM  M=  6 ,A)  , where  Q^,  and  A are  defined  as  above. 

The  transition  function  •+■  2 ^ maps  10-tuples  of  states 

into  sets  of  states.  That  is,  each  cell  has  nine  neighbors  — 
four  son  cells  in  a two-by-two  block  in  the  layer  below,  four 
nearest  neighbors  in  its  own  layer,  and  one  father  cell  in 
the  layer  above.  More  precisely,  let  M be  the  (i,j,k)th  cell 
in  a PA,  where  M is  in  the  ith  row  and  jth  column  of  retina 
cells  in  the  kth  layer  from  the  bottom.  Then  M's  father  is 
cell  ( f i/21 , f j/21 ,k+l) , where  fxl  denotes  the  smallest  integer 
greater  than  or  equal  to  x.  M's  brother  cells  are  (i-l,j,k), 
(i,j+l,k),  (i+l,j,k)  and  (i,j-l,k),  and  its  son  cells  are 
(2i,2j,k-l),  (2i,2j-l,k-l) , (2i-l , 2 j ,k-l)  and  (2i-l , 2 j-1 ,k-l) . 
The  figure  in  Section  1 illustrates  this  neighborhood. 

A configuration  of  a PA  is  an  assignment  of  states  from 
Qjj  to  each  cell  in  the  PA.  An  input  configuration  for  a 
pattern  of  size  2*^  by  2*^  is  an  n+3  high  stack  of  CA's.  The 
bottommost  layer  is  a size  2*^^^  by  2^^^  CA,  where  every  cell's 


initial  state  is  the  boundary  state  #.  The  next  layer  up  is 
called  layer  0 and  is  a size  2’^+l  by  2*^+1  CA,  whose  border 
cells  have  initial  state  #,  while  the  other  cells  define  the 
base  array  where  the  input  pattern  is  "stored"  — that  is, 
their  initial  states  are  chosen  from  and  define  the  input 
image  to  be  recognized.  Layer  1 is  a size  2'^~^+l  by  2'^~^+l  CA, 
with  retina  cells  initialized  to  the  quiescent  state  b£Q^. 
Layers  2 through  n are  similar  to  layer  1 except  each  layer's 
retina  has  one-quarter  the  number  of  cells  in  the  layer  below 
it.  Above  layer  n is  a single  cell  with  initial  state  # 
which  is  the  father  cell  of  the  apex  cell  in  layer  n. 

The  height  of  a PA  is  the  length  of  the  shortest  path 
from  the  root  to  a cell  in  the  base  array.  Thus  a PA  with 
base  size  2^^  by  2^  has  height  n and  4^+4””^+...+!  = (4*'^^-l)/3 
retina  cells.  An  input  image  is  a 2*^  by  2^  array  of  states 
from  for  some  niO,  which  defines  the  initial  states  of 
the  retina  cells  in  the  base  array  (layer  0) . A step  of 
computation  consists  of  a single  state  transformation  performed 
simultaneously  at  each  cell.  An  input  configuration  is  accepted 
by  a PA  if  the  root  enters  an  accept  state  after  a finite 
number  of  steps.  The  language  accepted  by  a PA  is  the  set  of 
all  input  images  accepted  by  it.  (As  with  CA's,  the  language 
is  defined  for  all  legal  input  sizes  and  is  not  associated 
with  a particular  instance  of  a PA  which  has  a specific  height 
and  recognizes  only  images  of  a specific  size.) 


A bottom-up  pyramid  cellular  acceptor  (BPA)  is  a PA  whose 

5 

state  transition  function  is  modified  to  be  6;Qj,  2 In 

N 

this  simplification  of  PA's,  the  next  state  of  a cell  depends 
only  on  the  current  states  of  that  cell  and  its  four  sons. 

Hence  information  can  only  move  up  the  pyramid.  All  other 
aspects  of  the  PA  hold  for  a BPA. 

Finally,  we  define  a second  simplification  of  PA's  which 
gives  another  alternative  neighborhood  definition  for  restrict- 
ing information  transmission  through  a cellular  pyramid.  An 

up-down  pyramid  cellular  acceptor  (UDPA)  is  a BPA  whose  state 

6 

transition  function  is  defined  as  6;Q„  ■*  2 . Here  the  next 

State  of  a cell  also  depends  on  the  state  of  its  father,  so 
that  state  information  can  move  up  and  down,  but  not  sidewise, 
through  the  pyramid. 

The  purpose  of  investigating  these  simplifications  of  the 
original  model  of  computation  is  not  because  of  any  hardware 
considerations  for  simplifying  the  interconnection  links. 

Rather,  it  is  intended  as  an  aid  in  studying  the  tradeoffs 
between  neighborhood  size,  time  and  space  (state  set  cardinality) 
bounds.  By  comparing  variations  of  the  model  we  gain  insight 
into  how  computing  power  is  affected  by  incrementally  adding 
more  channels  for  information  flow. 


3. 


BPA  Algorithms 


This  section  considers  several  basic  tasks  for  BPA's, 
including  detecting  and  counting  arbitrary  local  patterns, 
and  recognizing  palindromes,  rectangles  and  squares. 

We  assume  from  now  on  that  the  size  of  the  pyramid  base 
is  2^  by  2*^,  so  that  the  height  is  n.  The  layers  will  be 
numbered  0,1,..., n starting  from  the  base.  To  simplify  the 
exposition,  BPA's  will  not  be  defined  in  terms  of  states  and 
transition  functions;  rather,  we  will  specify  algorithms  in 
terms  of  transmitting  and  receiving  information  between  son 


and  father  cells. 


3 . 1 Local  Property  Detection 


In  this  section  we  consider  the  problem  of  detecting 

the  presence  or  absence  of  a d by  d pattern*  in  the  input  of 

a BPA.  The  difficulty  of  this  problem  is  that  such  a pattern 

can  be  badly  positioned  with  respect  to  the  pyramid  so  that 

only  cells  at  an  unbounded  height  above  the  base  can  see  all 

of  it.  Furthermore,  since  a d by  d block  of  the  base  must  be 

matched  against  the  given  pattern  at  the  lowest  common  ancestor 
2 

of  its  d base  cells  (because  this  cell  may  be  the  apex  cell 
of  the  BPA) , the  number  of  matchings  to  be  performed  by  a 
cell  increases  exponentially  with  its  height.  That  is,  a cell 
C in  layer  k must  check  all  those  d by  d blocks  containing 
base  cells  which  extend  across  the  "cracks"  between  c's  sons' 
bases.  The  figure  below  illustrates  the  base  cells  of  C which 
must  be  considered  in  the  matching  process  at  this  cell.  It 


*Since  any  arbitrary  finite  pattern  can  be  padded  with  "don't 
care"  symbols  to  obtain  an  equivalent  square  pattern,  there 
is  no  loss  in  generality  in  making  this  restriction. 


r 


is  easily  verified  that  there  are  2 (d-1) 2^  - 3 (d-1) ^ distinct 
d by  d blocks  to  be  checked  by  C. 

A one-dimensional  solution  to  this  problem  was 
given  in  [1] , but  that  solution  does  not  generalize  to  two 
dimensions.  In  this  section  we  present  a two-dimensional 
solution. 

Let  denote  the  first  (leftmost)  d-1  symbols  in 
cell  C's  ith  row,  and  let  denote  the  last  (rightmost)  d-1 
symbols  in  C ' s ith  row.  First,  construct  the  BPA  sc  that 

Jr 

cells  in  layer  k count  modulo  2 and  output  alternating 
k-1 

length-2  sequences  of  O's  and  I's  [1]. 

Suppose  cell  C at  height  6= flog  dl  repeatedly 
receives  a copy  of  the  entire  base  beneath  it  and 

a)  decides  whether  or  not  the  pattern  is  entirely 
contained  with  it. 

and  b)  modulo  2°  starting  at  time  6 outputs  and  C| 

for  1 s.  i ^ 2^ . 

Then  a cell  C in  layer  6+1  can  use  its  modulo  2*^^^  counters 
and  the  outputs  of  its  four  sons  to  compute  its  own  and 
C values  (modulo  2'^^^),  since 

r [C  (NW)  ] . , 1 a.  i_i  2^ 

UC(SW)1^_26  , 2^  < i 


t 


and 


[C(NE)  ][  , 1 ii  i2*^ 

C!  = - 
1 

, [C(SE)  ]]^_2<S  , 2*^  < i s.2'^'^^ 


where  C (MW)  is  C's  northwest  son,  etc. 

Simultaneously,  C can  decide  whether  or  not  the 
pattern  appears  anywhere  within  its  base  as  follows.  Since 
we  need  only  check  d-by-d  blocks  which  are  not  entirely  con- 
tained in  one  of  C's  son's  bases,  we  need  only  reconstruct  the 
2'^^^  by  2(d-l)  vertical  band  centered  on  the  middle  column  of 
C's  base  in  order  match  across  the  vertical  crack.  In  fact, 
only  a d by  2(d-l)  running  window  of  the  band  needs  to  be 
stored  by  cell  C at  any  time.  Using  its  modulo  2*^^^  counter 
and  its  four  sons'  outputs,  C can  scan  its  vertical  band 
since  the  ith  row  of  the  band  is  computable  at  time  6 + i 
(modulo  2*^^^)  by 

[<^(NW)]J1  [C(NE)]^  if  la.  is.  2^ 

or  [C(SW)  [C{SE)  if  2^^  < i a.  2^"^^ 

where  ||  denotes  concatenation.  After  time  d+6  C can  begin 
checking  whether  or  not  the  entire  pattern  is  present  anywhere 
in  its  d by  2{d-l)  window. 


In  the  same  way  C can  simultaneously  scan  the 
2(d-l)  by  2'^^^  horizontal  band  centered  on  the  middle  row  of 
C's  base  in  order  to  match  across  the  horizontal  crack.  (This 
means  that  the  cells  at  height  6 will  also  have  to  output  the 
top  and  bottom  d-1  symbols  in  the  ith  column,  1 i & 2 , 
modulo  2'^.)  If  at  any  time  the  pattern  is  found  in  either 
band,  a success  signal  is  immediately  sent  up  the  BPA  to  the 


apex  cell, 


Similarly,  each  cell  in  layer  k > 6+1  behaves  like 


the  cells  in  layer  6 + 1,  sequentially  scanning  their  length 

2 vertical  and  horizontal  bands  starting  at  time  k.  Hence, 

if  the  pattern  is  detectable  in  layer  k then  the  apex  cell 

K.  k 

will  know  this  at  a time  no  later  than  k + 2 + (n-k)  = 2 +n. 

In  the  worst  case  the  pattern  is  detectable  only  by  the  apex 

cell,  and  time  2'^t  (d)  + nt,  (d)  is  required,  where  t (d)  is 

the  time  needed  to  update  two  d by  2 (d-1)  windows  and  search 

for  the  d by  d pattern  within  them,  and  tj^(d)  is  the  time 

needed  to  copy  four  length  d vectors. 

If  we  assume  that  the  d by  d pattern  occurs  at 

least  once  and  is  equally  likely  to  appear  anywhere  in  the 

2*^  by  2^  base  and  d < < n,  the  0 (diameter)  worst  case  bound 

given  above  reduces  to  O(log  dieunater)  time  in  the  average 

case.  That  is,  there  are  cells  in  layer  k,  and  each 

k 2 

checks  2 (d-1) 2 - 3 (d-1)  d by  d blocks  for  the  presence  of 

the  pattern.  Thus  in  layer  k (k  > log  d)  a total  of 

2^"^  [(2  (d-1)  2*^  - 3{d-l)  ^)/2^’^l  d by  d blocks  are  inspected. 


It  is  easily  verified  that  after  processing  the  (6+2)nd 
layer  over  half  of  the  2^^  blocks  have  already  been  checked. 
Since  the  apex  cell  will  know  whether  or  not  the  pattern  is 
detected  in  any  of  these  layers  by  time  4d  t (d) 

Si 

the  0 (log  diameter)  average  time  bound  results. 


+ n tj^(d)  , 


Local  Property  Countin' 

Counting  the  number  of  occurrences  of  some  particular 
symbol  z in  a BPA's  base  is  an  immediate  generalization  of  the 
one-dimensional  counting  algorithm  [1] , which  required  0 (log 
diameter)  time.  The  problem  of  counting  arbitrary  local 
patterns  is  more  difficult  because  cells  at  all  levels  must 
detect  and  count  instances  of  the  pattern  and  then  add  their 
counts  in  with  all  other  cells'  counts. 

First,  we  show  that  a BTA  can  count  the  number  of 
occurrences  of  an  arbitrary  local  property  in  2 log  diameter 
time  steps.  That  is,  the  apex  cell  will  output,  at  time  steps 
n through  2n,  the  n+1  bits  in  the  binary  representation  of 
this  number,  least  significant  bit  first. 

Suppose  that  the  desired  pattern  has  length  d=2  . 

Each  cell  in  layer  6 receives  at  time  6 a copy  of  the  entire 
portion  B of  the  base  below  it  and  outputs,  at  the  end  of 
step  6, 

(1)  a 1 if  B matches  the  pattern,  and  a 0 otherwise; 

(2)  the  initial  d-1  length  segment  of  B,  call  this 
string  L; 

and  (3)  the  terminal  d-1  length  segment  of  B,  call  this 
string  R. 

Now  consider  a cell  C in  the  (k+l)st  layer.  By  the  induction 
hypothesis,  it  receives  from  its  sons  at  the  ends  of  time 
steps  k,...,2k  the  number  of  instances  of  the  given  pattern 


and  the  values  of  L and  R in  their  base  segments.  Let 


be  the  values  for  C's  left  son,  and  L2'^2  values  for  C's 

right  son.  Then  and  output  by  C at  the 

end  of  time  step  k+1.  In  addition,  Rj^||  is  the  length 
2(d-l)  segment  centered  on  the  middle  of  C's  base.  At  most 
d-1  occurrences  of  the  pattern  can  extend  across  C's  mid-base 
cell  and  all  the  information  for  checking  these  d-1  positions 
is  contained  in  the  string  R^||  L2 . Hence  at  time  k+1  C can 
compute  and  store  the  number  of  times,  S,  that  the  pattern  is 
found  crossing  C's  midpoint  and  then  add  the  least  signifi- 
cant bit  of  S to  the  sum  of  the  least  significant  bits  of  its 
sons'  counts,  which  has  just  been  computed.  At  step  k+2  C 
computes  its  next  least  significant  bit  by  adding  the  next 
bits  from  its  sons  and  the  next  bit  from  S.  This  process  is 
repeated  at  steps  k+3,...,k+6,  since  S is  a 6-bit  value.  For 
the  remaining  steps  k+6+1 , . . . , 2k+2  C just  sums  the  next  least 
significant  bits  from  its  sons.  Clearly  C's  outputs  at  steps 
k+l,...,2k+2  = 2 (k+1)  are  just  the  bits  of  the  sum  of  its 
sons'  counts  and  its  own  count,  i.e.,  the  number  of  times  the 
given  pattern  occurs  in  C's  base  segment.  The  total  time 
required  by  the  algorithm  is  2n  t +nt.  (d),  where  t is  the 
time  needed  to  perform  bit  addition  and  is  the  time 

required  to  count  the  nxamber  of  matches  of  a length  d pattern 
in  a length  2 (d-1)  string. 

This  algorithm  does  not  generalize  to  two  dimensions 
because,  as  we  noted  in  Section  3.1,  the  number  of  match  tests 
to  be  performed  at  a cell  is  now  unbounded.  Hence  a cell 


cannot  store  this  number  in  order  to  bit-serially  add  this 
count  in  with  the  counts  of  its  sons. 

By  allowing  sidewise  transmission  in  the  base,  counting 
arbitrary  local  patterns  reduces  to  the  problem  of  counting 
occurrences  of  a particular  symbol, since  the  base  array  can 
detect  and  mark  all  occurrences  in  2(d-l)  steps.  Similarly, 

other  algorithms  which  require  information  from  cells  only  a 
bounded  distance  apart  can  profit  from  sidewise  transmission 
in  the  base,  which  requires  only  a bounded  amount  of  time. 

For  example,  we  can  compute  the  Euler  number  of  a binary  image 
in  O(log  diameter)  time  on  such  an  extended  BPA  because  it 
is  a property  that  is  computable  from  measurements  taken  on 
2 by  2 windows  [2] . On  the  other  hand,  counting  the  number 
of  connected  components  in  a binary  image  using  a parallel 
shrinking  process  [3,4]  in  the  base  array  to  mark  components, 
and  then  using  the  upper  layers  of  the  pyramid  to  count  these 
marks,  provides  no  advantage,  since  the  shrinking  alone 
requires  up  to  O (diameter)  time. 


3.3  Palindromes 


In  two  dimensions  the  palindromes  may  be  defined 
in  alternative  ways  depending  on  the  axis  of  symmetry  chosen. 
Vertically  symmetric  palindromes  are  defined  as  those  arrays 
which  are  symmetrical  about  their  central  column.  That  is, 
if  the  input  symbols  of  a 2*^  by  2^  array  A are  indexed  in 
row-major  order  and  A (1 , 1)  =A (1 , 2*^)  , A (1 , 2)  =A (1 , 2*^-1)  , . . . , 
A(1,2”‘^-1)=A(1,2'^"^)  , . . . ,A(2'^,2’^"^-1)=A(2^,2"~^)  then  A is 
a palindrome. 

Before  presenting  the  algorithm  we  first  show  how 

ck+d 

a cell  in  layer  k of  a BPA  can  be  made  to  count  modulo  ab 
for  arbitrary  natural  numbers  a,b,c,d.  Since  ab^^^^=a  (b'^)  (b^)  ^ 
we  need  only  show  how  a cell  in  layer  k can  have  a modulo  ef 
counter  for  constants  e and  f.  We  have  shown  previously  [1] 
how  such  a counter  is  built  when  e=l  and  f=2.  In  the  general 
case,  base  cells  are  defined  to  output  I's  every  e time  steps, 
and  non-base  cells  output  I's  every  fth  time  that  their  sons' 

counters  output  I's.  Thus  layer  0 (base)  cells  are  modulo 

0 k 

ef  =e  counters.  If  cells  in  layer  k are  counting  modulo  ef  , 

then  readily  layer  k+1  cells'  f-step  delay  makes  these  cells 

f(ef^)=ef^^^  counters. 

Recognizing  palindromes  with  a BPA,  as  in  the  one- 
dimensional case,  is  based  on  performing  a fixed  sequential 
scan  of  its  base.  This  is  implemented  by  having  cells  in  the 

Jr 

kth  layer  count  modulo  4 as  described  above.  The  counters 

k-1 

can  be  easily  modified  so  as  to  output  0 for  the  first  4 


I 


k— 1 

steps,  1 for  the  next  4 steps,  2 for  the  next,  and  3 for 
k-1 

the  last  4 steps.  If  a father's  counter  is  0 it  copies 

its  NW  son's  value,  if  it  is  1 it  copies  its  NE  son's  value, 

if  it  is  a 2 it  copies  its  SW  son's  value,  and  if  it  is  a 3 

it  copies  its  SE  son's  value.  Thus  a cell  in  layer  1 copies 

its  sons'  values  at  successive  steps  in  the  following  order: 

^ . A cell  in  layer  2 copies  the  values  of  its  NW  son  for 
2 3 

four  steps,  then  its  NE  son  for  four  steps,  etc.,  so  that  it 
repeatedly  copies  the  values  of  its  base  cells  in  the  following 
order: 

0 14  5 

2 3 6 7 

8 9 12  13 

10  11  14  15 

By  induction,  it  follows  that  a cell  in  layer  k copies  the 
values  of  its  base  cells  at  times  k,k+l, . . . ,k+4  -1;  and  this 

If 

process  repeats  modulo  4 . 

We  can  obtain  a mirror  image  of  the  above  scan  by 
simply  changing  the  order  of  scanning  one's  sons  to:  ^ 2 • 

It  is  readily  seen  that  running  both  scans  simultaneously 
means  the  apex  cell  copies  all  pairs  of  symbols  in  base  cells 
which  are  vertically  symmetric  about  the  central  column  of 
the  base.  If  all  point  by  point  comparisons  match,  then  the 
BPA  can  accept  its  input.  The  total  time  required  is  n+4'^-1. 
This  is  nearly  optimal,  since  it  can  be  shown  that  the  lower 
bound  for  recognition  of  palindromes  on  a BPA  is  4^/2  by 


1 


remarks  similar  to  those  in  [1] . Alternative  definitions  of 
two-dimensional  palindromes,  for  example  by  using  horizontal 
or  central  symmetry,  can  be  handled  by  essentially  the  same 
method.  A related  question  which  is  still  open  is  whether 
the  apex  cell  of  a BPA  can  copy  the  symbols  in  its  base  in 
raster  order. 


■f 


3.4  Rectangles  and  Squares 


Let  the  language  consist  of  the  class  of 

input  images  which  contain  a single  solid  rectangle  of  I's 
(with  sides  parallel  to  the  sides  of  the  base)  on  a back- 
ground of  O's.  Blum  and  Hewitt  have  shown  [5]  how  an  FSA 
on  a two-dimensional  tape  can  recognize  by  checking 

each  boundary  point's  local  slope.  Minsky  and  Papert  (2] 
define  a diameter-limited  perceptron  to  recognize  by 

counting  the  local  property  "corner."  That  is,  input  image 

® ^ ^RECT 

(the  number  of  occurrences  of  pattern  ^ 

♦ (the  number  of  occurrences  of  pattern  ^ 4 

where  the  patterns  also  include  all  90°  rotations  of  each. 

A BPA  is  not  well  suited  either  for  simulating  the 
moves  of  an  FSA  (see  Section  4)  or  for  counting  local  properties 
(Section  3.2).  [If  we  allowed  the  BPA  to  have  sidewise  trans- 
mission in  its  base,  then  the  Minsky  and  Papert  method  could 
be  implemented  on  a BPA  and  would  require  0 (log  diameter) 

time  to  be  recognized.]  We  now  show  how  a BPA  can  recognize 
L__^_  by  a third  method  which  compares  consecutive  row  cross 
sections  of  the  base. 


Each  row  of  an  input  image  contained  in  must 

IT  s t 

be  of  the  form  010.  we  shall  now  prove  that  each  cell  C 
in  layer  k can  output  at  time  steps  k to  2k  the  k+1  bits  in 
the  binary  representations  of  r,s,  and  t for  the  top  row  and 


t 
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u,v,  and  w for  the  bottom  row  in  C's  subbase.  Associated 
with  the  first  (least  significant)  bit  of  each  number  will 
be  a flag  indicating  whether  or  not  the  value  is  zero. 

Clearly  the  cells  in  layer  1 can  compute  this 
information.  Now  consider  a cell  C in  the  (k+l)st  layer. 

By  the  induction  hypothesis,  it  receives  from  its  sons  at  the 
ends  of  time  steps  k,...,2k  the  values  of  r,s,...,w  in  their 
base  segments.  Let  the  subscripts  1,...,4  denote  the  north- 
west, northeast,  southwest,  and  southeast  sons  of  C,  respectively. 
Then  C's  base  has  the  following  form 


row  1 


row  2 


row  2 +1 

-k+1 
row  2 


so  that  C computes  its  own  r,  s and  t values  as  follows: 

1)  if  s^=0  and  S2=0  then  r=rj^+r2»  s=0,  t=0 

2)  if  s^^O  and  S2=0  then  s=Sj^,  t=tj^+r2 

3)  if  s^=0  and  S2?^0  then  r=rj^+r2»  s=S2,  t=t2 

4)  if  Sj^?<0  and  S2^0  and  tj^=0  and  ^2"®  then  ^=^2^,  s=Sj^+S2, 
t=t2 

If  none  of  these  conditions  holds,  a failure  signal  is 
immediately  propagated  up  to  the  root.  The  addition  is 
performed  as  in  the  BTA  counting  algorithm  (Section  4.2  of 


(1])  where  C functions  as  a bit  serial  adder.  Similarly,  C 


computes  u,  v,  and  w from  u^ , f ,w^ . 

Simultaneously,  C checks  whether  the  two  middle 
rows  in  its  base,  as  specified  by 


X V z 

match,  are  of  the  form  0 l-'O  and  are  consistent  with  its 

top  and  bottom  rows.  Say  these  two  middle  rows  are  of  the 
Xi  z X y z 

form  010  and  010,  respectively.  Then,  if  y^^O 
and  y2^0,  we  must  be  considering  two  rows  intersecting  a 
potential  rectangle.  Consequently,  only  if  Xj^=X2,  y2^~y2' 
and  can  both  rows  be  part  of  a legal  rectangle.  If 

yj^=y2=0,  both  rows  aro  scanning  only  background;  if  only  one 
of  y^^  and  y^  is  nonzero,  we  are  seeing  either  the  top  or 
bottom  of  the  rectangle.  In  these  cases  there  is  no  matching 
to  be  performed,  only  a consistency  check  that  prevents  multi- 
ple rectangles  from  being  accepted.  Details  will  not  be  given. 

By  transitivity  and  the  induction  hypothesis  we 
now  know  that  all  nonzero  rows  in  C's  base  segment  are 
connected  and  have  the  same  description,  since  every  consecu- 
tive pair  of  rows  within  C's  subpyramid  has  been  matched. 

If  a cell  ever  detects  a completed  rectangle,  a success  signal 
is  propagated  to  the  apex  cell  which,  if  it  receives  exactly 
one  such  signal,  accepts  its  input.  The  total  time  required 
by  the  algorithm  is  2nt  where  t is  time  needed  to  perform 


one  step  of  the  serial  addition  process  and  the  "sewing" 
together  of  quadrants  by  matching  boundary  rows. 

We  now  indicate  how  a BPA  can  detect  the  fact  that 
the  rectangle  of  I's  in  its  base  is  a square.  Define  L to 
be  a set  of  inputs  containing  a single  solid  square  of  I's 
on  a background  of  O's.  can  be  recognized  by  an  FSA  which 

first  checks  that  the  input  contains  a single  rectangle  and 
then  scans  from  this  rectangle's  upper-left  corner  at  45° 
toward  its  opposite  corner.  If  the  opposite  corner  is  reached, 
the  rectangle  is  a square,  otherwise  it  is  not  [5] . L 
cannot  be  recognized  by  a diameter-limited  perceptron  which 
measures  local  properties.  However,  there  does  exist  an 
order-limited  perceptron  for  recognizing  squares  which  works 
by  stratification  — a process  of  sequentially  enumerating 
and  testing  all  possible  starting  positions  and  side  lengths 
for  a square  in  the  image  [2] . Again,  neither  of  these  methods 
is  adaptable  for  use  on  a BPA. 

A BPA  can  recognize  L in  0 (log  diauneter)  time  by 
a modification  of  its  own  rectangle  algorithm.  Each  cell  C 
computes  cross  sections  of  its  base  segment's  boundary  rows 
and  columns.  In  addition,  if  C's  base  is  nonzero  C computes 
a row  and  column  cross  section  of  the  object  contained  in 
its  base.  This  is  done  as  follows:  if  C's  top  and  bottom 
rows  have  just  been  computed  to  be  all  zeros  and  C's  middle 
two  rows  are  nonzero,  then  C's  base  bounds  the  height  of  the 
object.  If  the  object  is  a rectangle,  C outputs  one  of  its 
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middle  rows,  which  is  a row  cross  section  of  the  rectangle. 

All  cells  above  C just  copy  this  cross  section  while  computing 
their  own  boundary  cross  sections.  Similarly,  cells  save  a 
column  cross  section  if  a rectangle's  width  is  bounded  by 
their  base  segment. 

When  a cell  detects  a completed  rectangle  it 
simultaneously  compares  the  I's  counts  in  the  horizontal  and 
vertical  cross  sections  through  the  rectangle.  If  they  are 
equal,  a success  signal  is  transmitted  to  the  apex  cell, 
otherwise  a reject  signal  is  sent.  The  apex  cell  accepts  its 
input  image  if  it  receives  exactly  one  accept  signal  and  no 
reject  signals  at  time  step  2n.  The  total  time  required  by 
this  algorithm  is  2n  t,  , where  t is  the  time  necessary  to 
update  six  cross  sections  and  "mend"  together  adjacent  quadrants 
by  matching  their  common  boundaries. 

Counting  the  number  of  rectangles  in  an  input  image 
by  a BPA  involves  the  same  problem  encountered  in  local  property 
counting  (Section  3.2).  That  is,  there  may  be  an  unbounded 
number  of  rectangles  which  are  detectable  only  by  the  apex  cell. 

The  general  technique  of  comparing  consecutive  row 
cross  sections  used  in  this  subsection  is  adaptable  for 
recognizing  other  fixed-orientation  polygons  whose  boundary 
slopes  are  locally  testable.  For  example,  straight  line  seg- 
ments with  slopes  which  are  multiples  of  45°  (with  respect  to 
the  bottom  edge  of  the  image)  are  digitally  realized  as 
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/ m 


T DOCUMENTATION  PAGE 


H’l-  AD  IN',  I NIK  TlriN^ 


X XX 

XXX  X X X 

X XX 

0°  45°  90°  135° 

Therefore,  a BPA  can  recognize  45°  right  triangles  with  one 
side  parallel  to  a side  of  the  base,  or  diamonds  (squares 
rotated  45°) , by  comparing  triples  of  consecutive  rows.  (Com- 
puting the  description  of  the  top  and  bottom  two  rows  of  each 
2 by  2 block  is  sufficient  for  this  test.)  Similarly, 
straight  lines  at  other  fixed  orientations  relative  to  the 
base  are  characterized  by  periodic  digital  segments  of  bounded 
length.  A BPA  can  be  constructed  to  detect  these  segments' 
slopes  by  examining  a bounded  number  of  cross  sections  above 
and  below  each  input  row. 

It  is  an  open  question  whether  a BPA  can  accept 
the  set  of  binary  images  whose  I's  form  a convex  subset. 
However,  several  necessary  conditions  for  an  image  to  be 
convex  can  be  tested  by  a BPA.  A binary  image  is  row  convex 
if  it  is  connected  and  each  row  contains  at  most  one  connected 
component.  This  is  a property  that  can  be  accepted  by  a BPA 
since  it  only  requires  verifying  that  there  is  at  most  one 
run  of  I's  on  each  row,  and  that  the  runs  of  I's  on  any  two 
adjacent  rows  overlap.  Similarly,  a BPA  can  recognize  column 
or  diagonal  convexity  . 


4 . The  Power  of  PA's 

In  this  section  we  investigate  the  power  of  PA's,  BPA's 
and  UDPA's  by  comparing  their  language  recognition  capabilities 
with  those  of  CA's  and  FSA's.  We  show  that  PA's,  UDPA's  and 
CA's  are  all  equivalent.  We  leave  open  the  question  of 
whether  or  not  BPA's  are  stronger  than  FSA's  on  a two- 
dimensional  tape. 


4 . 1 PA*s  and  CA*s 

In  [1]  we  showed  how  a one-dimensional  CA  can 
simulate  a one-dimensional  PA,  each  simultation  step  requiring 
O (diameter)  time.  This  result  i/nmediately  generalizes  to  two 
dimensions.  Clearly  PA's  can  simulate  CA's  since  only  cells 
in  the  base  array  have  boundary  cells  as  sons.  Thus  all  other 
cells  can  remain  in  the  quiescent  state  while  the  base  array 
copies  the  transitions  of  the  CA  using  only  state  information 
from  brother  cells.  If  the  CA's  upper-left  corner  cell  ever 
enters  an  accept  state,  the  father  of  the  upper-left  corner 
cell  in  the  PA's  base  array  can  detect  this  and  begin  propa- 
gating an  accept  signal  to  the  root.  Therefore,  PA's  are 
equivalent  to  CA's. 
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4.2  BPA's  and  FSA's 

We  have  previously  shown  that  BTA's  are  more  powerful 
than  FSA's,  and  can  simulate  them  in  Odog  diameter)  time . 

That  result  was  aided  by  the  fact  that  two-way  nondeterministic 
FSA's  are  no  more  powerful  than  one-way  deterministic  ones. 

In  two  dimensions,  however,  it  is  know  that  placing  restrictions 
on  the  allowable  directions  of  motion  does  alter  the  power  of 
FSA's.  (Henceforth,  FSA  will  mean  2-D  FSA.) 

An  FSA  on  a two-dimensional  input  tape  is  a 5-tuple 
M=  (Qj^,Q^,  (5 , A,qQ)  where,  as  in  the  one-dimensional  definition, 
is  a nonempty,  finite  set  of  states,  Qip  ^ Qj^  is  the  set  of 
input  states  (tape  symbols),  A c is  the  set  of  accept  states, 
and  qQ  Q is  the  start  state.  The  state  transition  function 
permits  four  directions  of  movement  — up,  down,  left  and  right. 

That  is,  6:Qj^  X X {U,D,L,R}  in  the  deterministic  case; 

(Q  x{U,D,L,R}) 

<5 : Qjg  X -*■  2 in  the  nondeterministic  case. 

To  begin  with  we  will  consider  a very  restricted 
type  of  FSA,  namely  one  that  can  only  do  a fixed,  "one-way" 
scan  of  its  tape.  In  particular,  we  show  how  a BPA  can  simu- 
late an  FSA  M that  can  only  do  a raster  (row-major)  scan  of 
its  input.  Tapes  which  are  not  square  and  whose  sizes  are  not 
powers  of  2 can  be  padded  with  a special  tape  symbol  3 at  the 
right  and  bottom  and  the  transition  function  altered  so  that 
M moves  over  }J's  without  changing  states.  We  assume  that  M 
is  deterministic,  since  clearly  one-way  nondeterministic  FSA's 
are  no  more  powerful,  by  the  same  arguments  used  to  show 


equivalence  in  the  one-dimensional  case  [6] . 

Let  the  state  set  of  M be  S={Sj^ , . . • , s^}  . Each 
cell  in  layer  k can  be  made  to  count  modulo  2 [1]  and  this 

counter  can  be  readily  modified  to  act  as  follows:  for  the 

k-1  k-1 

first  2 steps  it  outputs  0,  for  the  next  2 steps  it 

outputs  1,  and  so  on,  changing  from  0 to  1 or  1 to  0 every 

2^  ^ steps. 

At  time  step  1 each  base  cell  in  the  BPA  constructs 
a state  transition  vector  of  length  m based  on  the  input  value 
at  the  cell.  Each  non-base  cell  C behaves  as  follows:  when 
C's  counter  is  in  state  0 (1)  it  constructs  the  composition  of 
its  upper-left  (lower-left)  and  upper-right  (lower-right)  sons' 
transition  vectors.  Thus  a cell  in  layer  1 alternately  com- 
poses the  transition  vectors  of  its  upper  two  and  lower  two 
sons.  A cell  in  layer  2 composes  vectors  from  its  upper 
sons  at  the  first  two  time  steps  (i.e.,  the  first  two  rows  of 
its  base)  and  then  at  the  next  two  steps  computes  the  row  scan 
transition  vectors  for  its  third  and  fourth  rows.  By  induction, 
it  follows  that  a cell  in  layer  k computes  lef t-to-right  row 
scan  transition  vectors  for  the  2 rows  in  its  base  in  a top 

If 

to  bottom  sequence  at  times  k,k+l , . . . ,k+2  -1;  and  this  process 
repeats  (modulo  2 ) . 

At  the  same  time,  cell  C in  layer  k composes  the 
2 row  scan  transition  vectors  to  get  a complete  raster  scan 
transition  vector  for  C's  base.  This  is  accomplished  by  having 
C store  a second  cumulative  transition  vector.  When  C's 


f' 
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counter  switches  from  state  1 to  state  0 we  must  be  computing 
the  vector  for  C's  top  row.  Thus  at  this  step  we  initialize 
the  raster  scan  transition  vector  with  the  first  row's  vector. 
At  all  other  times,  we  just  compnase  the  current  row  p's  vector 
with  the  saved  raster  scan  vector  to  obtain  a new  raster  scan 
vector  describing  M's  movement  over  the  top  p rows  of  input. 

In  particular,  at  time  n+2”-l  the  apex  cell's  raster  scan 
transition  vector  describes  the  state  that  M ends  in  after 
scanning  the  base  in  raster  order  when  starting  in  any  of  its 
m states.  Thus  if  M's  initial  state  gives  rise  to  an  accepting 
state  in  this  vector,  the  BPA  can  accept,  otherwise  it  rejects 
its  input. 

If  the  time  required  to  look  up  a value  in  a table 

* 

of  length  k is  t(k),  then  t(lQ,pl)  is  the  time  necessary  for 
the  base  cells  to  initially  set  up  their  transition  vectors. 

An  additional  (n+2”-l)t(m)  time  is  then  required  before  the 
apex  cell  can  decide  membership  in  L (M)  . Thus  t(|Q^|)  -t- 
(n-f 2*^-1)  t (m)  or  O(diameter)  time  is  sufficient.  This  is 
faster  than  M itself,  which  requires  O (area)  time  to  scan  its 
input.  However,  a CA  can  also  simulate  a raster  scan  FSA  in 
0 (diameter)  time. 

Similarly,  we  can  simulate  other  fixed  scanning 
sequenoea  of  an  FSA,  e.g.,  snake-like  or  colvunn-major  indexed 
scans,  without  significantly  altering  the  construction  given 
above.  The  critical  knowledge  that  we  have  used  here,  which 
is  not  available  in  more  general  FSA's,  is  that  the  motion  of 


H is  fixed  in  advance  and  so  it  is  not  necessary  to  save  M's 
position  on  the  tape. 

We  now  return  to  the  general  problem:  how  does  a 
BPA  compare  with  an  FSA  which  can  move  to  any  of  its  four 
adjacent  neighbors  during  a transition?  The  problem  is  that 
M can  enter  or  leave  a given  block  of  base  cells  at  any  place 
along  its  boundary,  and  so  the  description  of  M's  behavior  on 
a block  grows  with  the  size  of  the  block.  Consequently,  we 
cannot  specify  M's  behavior  relative  to  an  unbounded  size  base 
segment  by  a state  transition  vector  of  bounded  length. 

The  search  for  an  alternative  method  is  worthwhile, 
however,  since  the  following  result  shows  that  a BPA  has 
sufficient  time  to  distinguish  between  all  possible  base 
segments.  Let  m be  the  number  of  states  in  M and  let  s ^m  be 

the  number  of  input  states.  Then  clearly  the  number  of  possible 

k k 2^^ 

2^  by  2^  input  configurations  is  D=s  . From  [1],  we  know 

If 

that  a cell  in  layer  k of  a BPA  (with  a 2 by  2 base  segment) 
can  have  a state  sequence  period  up  to  ml  (2*^)^°^  It  follows 

that  the  total  number  of  distinct  sequences  of  states  that  a 

m. ^ 

cell  in  layer  k can  have  is  bounded  from  above  by  P=m 

However,  when  k>>m  we  have  D <P,  implying  that  a BPA's  period- 

icity is  not  a limitation  on  distinguishing  between  all  possible 
2*^  by  2^  input  blocks. 

In  any  event,  the  sets  of  languages  accepted  by 
BPS's  and  FSA's  are  not  the  same.  The  vertically  symmetric 
palindromes  were  shown  in  Section  3.3  to  be  recognizable  by  a 


BPA.  Blum  and  Hewitt  have  proved  [5]  that  this  language 
cannot  be  accepted  by  an  FSA. 


4 . 3 UDPA's  and  PA's 

UDPA's  are  a simplification  of  PA's  in  which  the 
neighborhood  definition  eliminates  all  sidewise  connections 
to  brother  cells.  Though  slower  for  some  tasks,  we  show  that 
UDPA's  are  equivalent  to  PA's.  Clearly,  any  language  that 
can  be  recognized  by  a UDPA  can  be  recognized  by  a PA  that 
ignores  its  brother  links.  We  now  show  that  any  language 
recognized  by  a PA  can  also  be  recognized  by  a UDPA.  We 
prove  this  by  demonstrating  how  a UDPA  can  simulate  a CA. 

For  simplicity,  we  give  the  one-dimensional  proof;  the  generali- 
zation to  the  two-dimensions  is  immediate. 

Given  a CA  with  input  size  N,  let  n be  the  smallest 
integer  such  that  The  input  string  will  be  left- 

justified  in  the  base  of  the  UDPA  with  the  rightmost  2^^  - N 
cells  initialized  to  the  boundary  state  #. 

To  simulate  a single  step  of  the  CA,  each  UDPA  base 
cell  must  have  access  to  the  states  of  its  corresponding  CA 
cell's  two  brother  cells.  Every  two  cells  that  are  adjacent 
in  the  CA  have  a least  common  ancestor  (LCA)  in  the  UDPA. 

Furthermore,  that  UDPA  cell  which  is  the  LCA  for  an  adjacent 
pair  in  the  CA  cannot  be  the  LCA  for  any  other  pair.  Therefore, 
every  two  adjacent  cells  in  the  base  of  a UDPA  can  simul- 
taneously switch  state  information  as  follows; 

At  time  step  1,  cells  in  layer  1 copy  the  pair  of 
states  in  their  base  segments.  At  time  step  2,  left  (right) 
sons  in  layer  0 copy  the  right  (left)  state  stored  in  their 

f 


father's  cell.  At  the  same  time  cells  in  layer  2 store  the 
ordered  quadruple  of  states  in  their  base  segments  by  copying 
the  state  pairs  stored  in  layer  1 cells.  Thus  at  the  end  of 
step  1 N/2  LCA's  have  been  found,  and  at  the  end  of  step  2 
each  base  cell  knows  the  state  of  one  of  its  brothers.  At 


step  3,  left  (right)  sons  in  layer  1 copy  the  third  (second) 
state  stored  in  their  father's  quadruple  of  states  and  mark 
it  L(R).  Also,  cells  in  layer  3 compute  their  own  ordered 
quadruple  of  states  by  copying  the  first  and  fourth  members 
from  their  two  son's  quadruples.  Thus  at  the  end  of  step  2 
N/4  more  LCA's  are  found  in  layer  2.  After  step  3 the  states 
of  these  N/4  adjacent  cell  pairs  have  been  swapped  by  the  sons 
of  their  LCA's.  At  subsequent  steps,  any  right  (left)  son  of 
a father  with  state  marked  L(R)  copies  its  father's  state. 

In  this  way  the  state  information  propagates  back  down  the 
UDPA  to  the  proper  brother  cell. 

Similarly,  state  quadruples  are  computed  higher 
and  higher  in  the  UDPA  until  every  state  pair's  LCA  is  found. 
At  time  n the  root  is  left  with  a pair  of  states  corresponding 
to  the  states  of  the  leftmost  and  rightmost  cells  in  the  base. 
Consequently,  their  brothers'  states  are  the  boundary  state  #, 
and  this  information  can  be  returned  to  the  proper  base  cell 
in  the  same  manner.  Thus  simulating  one  step  of  a CA  takes 
2n  or  O(log  diameter)  time  for  the  UDPA. 


5. 


Conclusion 


We  have  extended  our  study  of  pyramid  cellular  acceptors 
in  this  report,  treating  primarily  the  acceptance  of  two- 
dimensional  properties  by  bottom-up  pyramid  acceptors.  In 
particular,  we  have  described  how  BPA's  can  do  local  property 
detection,  but  not  counting.  Since  PA's  can  count  local 
properties,  the  advantage  of  sidewise  connections  (at  least 
in  the  base)  has  been  established.  In  addition,  we  showed 
how  a BPA  can  accept  two-dimensional  palindromes,  a language 
which  is  not  recognizable  by  any  two-dimensional  FSA.  While 
PA's  are  shown  to  accept  precisely  the  class  of  languages 
accepted  by  CA's,  it  is  left  open  whether  BPA's  can  accept 
the  two-dimensional  finite  state  languages. 

Table  I summarizes  some  of  the  BPA  recognition  results 
and  open  questions  to  date  for  one  and  two-dimensional 
languages . 
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Recognition  Times  for  BPA's  and  BTA' 


Here  d is  the  string  length  or  array  diameter.  For  CA' 
the  corresponding  recognition  times  are  all  0(d). 
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